Magnetism and superconductivity in the t—t'—J model 
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We present a systematic study of the phase diagram of the t—t'—J model by using the Green's 
function Monte Carlo (GFMC) technique, implemented within the fixed-node (FN) approximation 
and a wave function that contains both antiferromagnetic and d-wave pairing. This enables us 
to study the interplay between these two kinds of order and compare the GFMC results with 
the ones obtained by the simple variational approach. By using a generalization of the forward- 
walking technique, we are able to calculate true FN ground-state expectation values of the pair- 
pair correlation functions. In the case of t' = 0, there is a large region with a coexistence of 
superconductivity and antiferromagnetism, that survives up to 8 C ~ 0.10 for J/t = 0.2 and S c ~ 0.13 
for J/t — 0.4. The presence of a finite t' /t < induces a strong suppression of both magnetic (with 
5 C < 0.03, for J/t — 0.2 and t 1 /t = —0.2) and pairing correlations. In particular, the latter ones are 
depressed both in the low-doping regime and around 8 ~ 0.25, where strong size effects are present. 

PACS numbers: 



I. INTRODUCTION 



After more than twenty years from the discovery 
of high-temperature superconductivity, a comprehensive 
description of the cuprate materials is still lacking. One 
of the main concern is about the origin of the electron 
pairing, namely if it is due to electron-phonon coupling, 
like in the standard theory by Bardeen, Cooper and 
Schrieffer (BCS), 1 or it can be explained by alternative 
mechanisms, based on the electronic interaction alone. 
From one side, though the isotope effect in cuprates (if 
any) is much smaller than the one observed in BCS su- 
perconductors, there are experiments suggesting a strong 
coupling between electrons and localized lattice vibra- 
tionsi^ On the other side, besides a clear experimental 
outcome showing unusual behaviors in both metallic and 
superconducting phases, there is an increasing theoretical 
evidence that purely electronic models can indeed sus- 
tain a robust pairing, possibly leading to a high critical 
temperature^ 5 -! 6 - Within the latter scenario, the minimal 
microscopic model to describe the low-energy physics has 
been proposed to be the Hubbard model or its strong- 
coupling limit, namely the t—J model, which includes an 
antiferromagnetic coupling between localized spins and 
a kinetic term for the hole motioni^ In this respect, 
Anderson proposed that electron pairing could naturally 
emerge from doping a Mott insulator, described by a res- 
onating valence bond (RVB) state, where the spins are 
coupled together forming a liquid of singlets^ Indeed, 
subsequent numerical calculations for the t—J model, 9 
showed that, though the corresponding Mott insulator 
(described by the Heisenberg model) has magnetic or- 
der, the RVB wave function with d-wave symmetry in 
the electron pairing can be stabilized in a huge region 
of doping close to the half-filled insulator. These calcu- 
lations have been improved by studying the accuracy of 
such a variational state, giving solid and convincing argu- 
ments for the existence of a superconducting phase in the 
t—J model. 4 However, other numerical techniques, like 



Density Matrix Renormalization Group (DMRG), pro- 
vided some evidence that charge inhomogeneities can oc- 
cur at particular filling concentrations ! 10 ! 11 These stripes 
are probably enhanced by the strong anisotropic bound- 
ary conditions used in this approach and can be also 
found by allowing anisotropies in the hopping and in 
the super-exchange coupling^ Coming back to the pro- 
jected RVB wave functions, it is worth mentioning that 
an approximate and simplified description of these states 
can be obtained by the renormalized mean-field theory 
(RMFT) , the so-called "plain vanilla" approach^ When 
this approach is applied to the t — J model, it is pos- 
sible to describe many unusual properties of the high- 
temperature superconductors and capture the most im- 
portant aspects of the cuprate phase diagram^ However, 
at present, most of the calculations have been done by 
neglecting antiferromagnetic correlations, that are defini- 
tively important at low doping. Within RMFT and most 
of the variational calculations, the magnetic correlations 
are omitted, implying a spin liquid (disordered) state in 
the insulating regime. Although antiferromagnetism can 
be easily introduced in both approaches, it is often not 
satisfactorily described, since the presence of an antifer- 
romagnetic order parameter in the fermionic determinant 
implies a wrong behavior of the spin properties at small 
momenta ) 15 ! 16 unless a spin Jastrow factor is used to de- 
scribe the corresponding spin-wave fluctuations. Indeed, 
it is now well known that the accurate correlated de- 
scription of an ordered state is obtained by applying a 
long-range spin Jastrow factor to a state with magnetic 
orde n 17 : 18 ' 19 The important point is that the Gaussian 
fluctuations induced by the Jastrow term must be orthog- 
onal to the direction of the order parameter, in order to 
reproduce correctly the low-energy spin-wave excitations. 
Moreover, by generalizing the variational wave function 
to consider Pfaffians instead of simple determinants ) 20 i 21 
it is possible to consider both electron pairing and mag- 
netic order, that are definitively important to determine 
the phase diagram of the t—J model in the low-doping 
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regime. 

The interplay between superconductivity and mag- 
netism is the subject of an intense investigation in the 
recent years. In most of the thermodynamic measure- 
ments these two kinds of order do not coexist, though 
elastic neutron scattering experiments for underdoped 
YBa 2 Cu30 x could suggest a possible coexistence, with 
a small staggered magnetization . 22 ' 23 ! 24 On the contrary, 
in the t— J model, there is an evidence in favor of a co- 
existence,— the antiferromagnetic order surviving up to a 
relatively large hole doping, i.e., 6 ~ 0.1 for J/t = 0.2. 21 
Therefore, the regime of magnetic order predicted by 
these calculations extend to much larger doping than 
the experimental results and also the robustness of the 
coexistence of superconductivity and antiferromagnetism 
seems to be inconsistent with the experimental outcome. 
Of course, disorder effects, which are expected to be im- 
portant especially in the underdoped region, would affect 
the general phase diagram^ However, without invoking 
disorder, one is also interested to understand if alterna- 
tive ingredients can modify the phase diagram of the sim- 
ple t—J model. For instance, band structure calculations 
support the presence of a sizable second-neighbor hop- 
ping if in cuprate materials, showing a possible connec- 
tion between the value of the highest critical temperature 
and the ratio t' /t£& Moreover, an experimental analysis 
suggests an influence of the value of t'/t on the pseudo- 
gap energy scaled From a theoretical point of view, the 
effect of t' is still not completely elucidated^^^Li 3 ^ 3 . 
though there are different calculations providing evidence 
that a finite t' could suppress superconductivity in the 
low-doping regime. On the other hand, recent Monte 
Carlo calculations suggest that the presence of if (as well 
as a third- neighbor hopping t") could induce an enhance- 
ment of pairing in optimal and overdoped regions. 31 ! 32 

In this paper, we want to examine the problem of the 
interplay between magnetism and superconductivity in 
the t—J model and its extension including a next-nearest- 
neighbor hopping t' by using improved variational and 
Green's function Monte Carlo (GFMC) techniques. In- 
deed, especially the latter approach has been demon- 
strated to be very efficient in projecting out a very accu- 
rate approximation of the exact ground state and, there- 
fore, can give useful insight into this important issue re- 
lated to high-temperature superconductivity. 

The paper is organized as follows: In Sec. |TT] we de- 
scribe the methods we used, in Sec. IIIII we show the nu- 
merical results for antiferromagnetism and superconduc- 
tivity, and in Sec. IIVI we draw our conclusions. 



II. MODEL AND METHOD 
A. Model and variational wave function 

We consider the t—t'—J model on a two-dimensional 
square lattice with L sites and periodic boundary condi- 



tions on both directions: 



- * J2 c U c ^ - X! c U ci '° + h - c - 

<i,J> ((k,l))<? 

where (...) indicates the nearest-neighbor sites, ((. . . )) 
the next-nearest-neighbor sites, cj a (ci )CT ) creates (de- 
stroys) an electron with spin a on the site i, 
Si = (S?,S?,Sf) is the spin operator, Sf = 
a' c l cr r er<T' c i,cr' > being r Q the Pauli matrices, and 
rii = cJ CT Cj jCT is the local density operator. In the fol- 
lowing, we set t = 1 and consider if = and t'/t = —0.2. 
Moreover, we consider two kinds of square clusters: Stan- 
dard clusters with L — I x I sites and 45° tilted lattices 
with L = 2 x I 2 sites. Besides translational symmetries, 
both of them have all reflection and rotational symme- 
tries. 

The variational wave function is defined by: 



\^VMC) = JsJdV N P G \$ 



MF 



(2) 



where Vq is the Gutzwiller projector that forbids double 
occupied sites, Vn is the projector onto the subspace 
with fixed number of N particles. Moreover, J s is a spin 
Jastrow factor 



(3) 



being Vij variational parameters, and finally Jd is a den- 
sity Jastrow factor 



J d = exp 



h3 



(4) 



being Uij other variational parameters. The above wave 
function can be efficiently sampled by standard varia- 
tional Monte Carlo, by employing a random walk of a 
configuration \x), defined by the electron positions and 
their spin components along the z quantization axis. In- 
deed, in this case, both Jastrow terms are very simple 
to compute, since they only represent classical weights 
acting on the configuration. 

As previously reported,— the main difference from pre- 
vious approaches is the presence of the spin Jastrow fac- 
tor and the choice of the mean-field state |$mf), defined 
as the ground state of the mean-field Hamiltonian 

i,j.<7 i.a 

+ a ^( c It c L + c 1t c u + h - c? > + Uaf i ( 5 ) 

(*,J> 

where we include both BCS pairing Ajj [with d- 
wave symmetry, i.e., for nearest-neighbor sites = 
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A(cosfca: — cosfcj,)] and staggered magnetic field in the 
x— y plane 



haf = 



(0) 



where Aaf is a variational parameter that, together with 
the chemical potential fi and the next-nearest-neighbor 
hopping of Eq. §E§ , can be determined by minimizing the 
variational energy of TL. Whenever both A and Aaf are 
finite, the projection (x\<$>mf) of the mean-field state on 
a given configuration \x) can be described in terms of 
Pfaffians, instead if A = or Aaf = it can be de- 
scribed by using determinants. Moreover, only in the 
case where the magnetic order parameter is in the x—y 
plane, the presence of the spin Jastrow factor ((3]) can 
introduce relevant fluctuations over the mean-field order 
parameter Aaf, leading to an accurate description of 
the spin properties. A detailed description of the wave 
function of Eq. (J2J) and its physical properties can be 
found in Ref. l2ll The variational parameters contained 
in the mean-field Hamiltonian (J5]) and in the Jastrow fac- 
tors ^ and ([3]) are calculated by using the optimization 
technique described in Ref. [H, that makes it possible to 
handle with a rather large number of variational param- 
eters. 



B. GFMC: beyond the Variational Monte Carlo 

The optimized variational wave function \^vmc) can 
be also used within the GFMC method to filter out an 
approximation of the ground state \^>q N ). Indeed, due 
to the presence of the fermionic sign problem, in order to 
have a stable numerical calculation, the GFMC must be 
implemented within the fixed-node (FN) approach, that 
imposes to ) to have the same nodal structure of the 
variational ansatzj 3 -^ Here, we recall the basic definitions 
of the standard FN method. A detailed description of 
this technique can be found in Ref. [2ll 

Starting from the original Hamiltonian TL, we define 
an effective model by 



TL 



eff 



TL + O. 



(J) 



The operator O is defined through its matrix elements 
and depends upon a given guiding function that is 
for instance the variational state itself, i.e., |\EVmc) : 



Ox' x — 



TLx r : x 

'V,»V,ai>0 



if Sx>,x = - ft, , > 
|* for x 1 = x, 



where ty x — { x \^}, \ x ) being an electron configuration 
with definite z-component of the spin. Notice that the 
above operator annihilates the guiding function, namely 
0\ty) = 0. Therefore, whenever the guiding function is 
close to the exact ground state of TL, the perturbation 
O is expected to be small and the effective Hamiltonian 
becomes very close to the original one. The most im- 
portant property of this effective Hamiltonian is that its 



ground state \^q N ) can be efficiently computed by us- 
ing GFMC.2&2! The distribution H x oc (x\^)(x\^ N ) is 
sampled by means of a statistical implementation of the 
power method: n oc lim IWOO G n U°, where n° is a start- 
ing distribution and G x >,x = *x'(A5 x ',a; - 'Heff,x',x)/'^x, 
is the so-called Green's function, 5 X >, X being the Kro- 
necker symbol. The statistical method is very efficient 
since all the matrix elements of G are non-negative 
and, therefore, G can represent a transition probability 
in configuration space, apart for a normalization factor 
b x = J2 X ' G X ',x- Since \^q N ) is an exact eigenstate of the 
effective Hamiltonian TL e ff, the corresponding ground- 
state energy can be evaluated efficiently by computing 



Ema 



V Mc\n eff \^ N ) 



(*VMCM N ) 



(8) 



namely the statistical average of the local energy e l (x) = 
{^vmc\T~L\x)/{ 1 &vmc\x) over the distribution n^. The 
quantity Em A < Evmc because, by the variational prin- 



ciple E 



ma 



< 



(¥|We//|*)/<*|¥) = E 



VMC- 



Moreover, 



Ema represents an upper bound of the expectation value 
Efn of TL over as it is shown in Ref. |35| or it 

can be simply derived by considering that the opera- 
tor O is semi-positive definite, namely all its eigenval- 
ues are non-negative. In the following, we will denote 
by FN the (variational) results obtained by using the 
GFMC method with fixed-node approximation, whereas 
the standard variational Monte Carlo results obtained by 
considering the wave function of Eq. @ will be denoted 
by VMC. 

Summarizing the FN approach is a more general and 
powerful variational method than the straightforward 
variational Monte Carlo. Within the FN method, the 
wave function 1^^), the ground state of 7i e ^ is known 
only statistically, and, just as in the variational approach, 
Efn depends explicitly on the variational parameters 
defining the guiding function \^}. This is due to the 
fact that TLeff depends upon |\&) through the operator 
O. The main advantage of the FN approach is that it 
provides the exact ground-state wave function for the 
undoped insulator (where the signs of the exact ground 
state are known), and therefore it is expected to be par- 
ticularly accurate in the important low-doping region. 
Moreover, the FN method is known to be very efficient in 
various cases: For instance, it has allowed to obtain a ba- 
sically exact description of the three-dimensional system 
of electrons interacting through the realistic Coulomb po- 
tential (in presence of a uniform positive background). 38 
Therefore, it represents a very powerful tool to describe 
correlated electronic systems. 



III. RESULTS 

A. Phase separation 

Before showing the results on magnetic and supercon- 
ducting properties, we briefly discuss the stability against 
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phase separation. In order to detect a possible phase sep- 
aration, it is very useful to follow the criterion given in 
Ref. [H and consider the energy per hole: 



e h {5) 



e(S) - e(0) 



(9) 



where e(S) is the energy per site at hole doping d and 
e(0) is its value at half filling. In practice, eh{5) repre- 
sents the chord joining the energy per site at half fill- 
ing and the one at doping 8. For a stable system, eh{8) 
must be a monotonically increasing function of 6, because 
the ground-state energy of a short-range Hamiltonian is 
a convex function of the doping. By performing exact 
energy calculations on finite clusters, the phase separa- 
tion instability is marked by a minimum of eh{8) at a 
given 6 C and by a flat behavior up to S c in the thermody- 
namic limit. For an approximate variational technique 
based on a spatially homogeneous ansatz, this flat be- 
havior is never reached in the phase separated region, so 
that eh{5) remains with a well defined minimum even for 
very large sizes; in this case, S c can be estimated with the 
Maxwell construction, provided the variational asnatz is 
accurate enough. In Ref. |2lj, we demonstrated the exis- 
tence of an homogeneous state for t' = and J/t < 0.7. 
As shown in Table HI the FN approximation, that is ex- 
act at zero doping^ provides a substantial lowering of 
the VMC energy, especially away from half filling and 
for a finite t' . This is a first indication that, for t'/t < 0, 
the simple variational approach could not be adequate to 
provide a reliable quantitative description of the ground- 
state properties. 

The FN results clearly indicate that the inclusion of 
a negative next-nearest-neighbor hopping contributes to 
further stabilize the homogeneous phase at finite doping, 
see Fig. [TJ This result is compatible with the outcome 
of recent calculations based on cluster dynamical mean- 
field theory on the Hubbard model, where a negative ra- 
tio t'/t enhances the stability of the homogeneous phase, 
whereas positive values of t' favor phase separation^ In 
this work, we do not want to address in much detail this 
issue and we will focus our attention on the more inter- 
esting magnetic and superconducting properties. 



B. Antiferromagnetic properties 

Here we present the results for the magnetic properties 
of the t— t'— J model and compare the FN approach with 
the VMC one, based upon the wave function ©. As 
already discussed in Ref. l2ll . the optimized wave func- 
tion @ breaks the SU(2) spin symmetry, because of the 
magnetic order parameter Aaf of Eq. and the spin 
Jastrow factor ([3]). It turns out that at half- filling and in 
the low-doping regime, the variational state ([2]) has an 
antiferromagnetic order in the x—y plane, whereas the 
spin-spin correlations in the z axis decay very rapidly. 
Therefore, in order to assess the magnetic order at the 
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FIG. 1: (Color online) Energy per hole eh{5), calculated by 
using the FN method, as a function of the doping 8 for L — 98 
and 162 and two values of the next-nearest-neighbor hopping 
t'/t = and -0.2. 



TABLE I: Variational (VMC) and fixed-node (FN) energies 
per site for J/t = 0.2 and t' = (third and fourth columns), 
and t'/t — —0.2 (fifth and sixth columns) for two clusters with 
L = 98 and 162 and different hole concentrations Nh — L — N 



L 


N h 


Evmc / L 


Efn/L 


Evmc / L 


Efn/L 


98 





-0.233879(1) 


-0.23432(1) 


-0.233879(1) 


-0.23432(1) 


98 


2 


-0.274144(5) 


-0.27752(1) 


-0.27290(1) 


-0.27808(1) 


98 


4 


-0.31429(1) 


-0.32053(1) 


-0.31189(1) 


-0.32123(1) 


98 


6 


-0.35482(1) 


-0.36328(1) 


-0.35132(1) 


-0.36405(1) 


98 


8 


-0.39550(1) 


-0.40563(2) 


-0.39028(1) 


-0.40575(1) 


98 


10 


-0.43581(1) 


-0.44728(2) 


-0.42814(1) 


-0.44561(1) 


162 





-0.233707(1) 


-0.23409(1) 


-0.233707(1) 


-0.23409(1) 


162 


2 


-0.258002(5) 


-0.26020(1) 


-0.257260(5) 


-0.26012(1) 


162 


4 


-0.282117(5) 


-0.28621(1) 


-0.28067(1) 


-0.28698(1) 


162 


6 


-0.306324(5) 


-0.31212(1) 


-0.30429(1) 


-0.31307(1) 


162 


8 


-0.33060(1) 


-0.33793(1) 


-0.32807(1) 


-0.33925(2) 


162 


10 


-0.35498(1) 


-0.36360(2) 


-0.35207(1) 


-0.36514(2) 


162 


12 


-0.37954(1) 


-0.38912(2) 


-0.37567(1) 


-0.39079(2) 


162 


14 


-0.40406(1) 


-0.41446(2) 


-0.39939(1) 


-0.41520(2) 


162 


16 


-0.42838(1) 


-0.43946(2) 


-0.42232(1) 


-0.43936(2) 



variational level, we have to consider the isotropic spin- 
spin correlations: 

la c\ (*VMcl S ' Srl^VMc) n n x 

(S>o-5s r ) = — . (10) 



VMCV&VMC 



The FN approach alleviates the anisotropy between the 
x—y plane and the z axis; in this case, we find that a 
rather accurate (and much less computational expensive) 
way to estimate the magnetic moment can be obtained 
from the z component of the spin-spin correlations: 



(*$ N \n N ) ' 



(ii) 



This quantity can be easily computed within the forward- 
walking technique^ because the operator SqS* is diago- 
nal in the basis of configurations used in the Monte Carlo 
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FIG. 2: (Color online) Magnetization obtained from the spin- 
spin correlations at the maximum distance calculated for the 
t— J model with J/t = 0.2 (upper panel) and J/t = 0.4 (lower 
panel). For the VMC calculations the error-bars are smaller 
than the symbol sizes. The VMC magnetization has been 
obtained from the isotropic correlations, whereas the FN one 
from the correlations along the z axis (see text). 



sampling. From the spin-spin correlations at the maxi- 
mum distance, it is possible to extract the value of the 
magnetization. In particular, the variational wave func- 
tion is not a singlet when the antiferromagnetic order sets 
in, and the magnetization has to be computed with the 
spin isotropic expression M — Hindoo y/ (So ■ S r ). On 
the other hand, the FN ground state is almost a perfect 
singlet for all the cases studied and the magnetization can 
be estimated more efficiently by M — linv^oo y/3{S^S^). 
The spin isotropy of the FN wave function can be explic- 
itly checked by computing the mixed-average of the total 
spin square 



(s 2 ) 



MA 



VMC 



s- 



I*™ 



VMC 



(12) 



that vanishes if |\&q ) is a perfect singlet, even when 
Y&vmc) has not a well-defined spin value. 

In Fig. [2] we report the results of the magnetization 
in the t—J model with J/t = 0.2 and 0.4. At finite 
doping, it is not possible to perform a precise size scal- 
ing extrapolation since it is very rare to obtain the same 
doping concentration for different cluster sizes. More- 
over, though the FN approach is able to recover an exact 
singlet state at half filling, (S 2 )ma increases by doping, 
reaching its maximum around 5 ~ 0.06, e.g., (S 2 ) ma ~ 1 
for 8 holes on 162 sites. This could explain why the FN 
results are a bit larger than the VMC ones for 6 ~ 0.06, 
especially for J/t = 0.2. Definitively, both the VMC and 
FN wave functions are almost spin singlets close to the 
transition point, because the mean-field order parameter 
Aaf goes to zero together with the parameters defining 
the spin Jastrow factor. Therefore, we are rather confi- 
dent in the estimation of the critical doping 8 C , where the 
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FIG. 3: (Color online) Static spin structure factor S(q) for 
L = 16 x 16 cluster and different hole concentrations for the 
t-J model with J/t = 0.2. T = (0,0), X = (tv,ty), and 
M — (7r,0). Inset: S(q) for the variational state (empty 
symbols) and for the FN approximation (full symbols). 
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FIG. 4: (Color online) Spin structure factor S(q) for the t—J 
with J/t = 0.2 and doping 5—1/8 and different clusters 
(L = 8 x 8, 12 x 12, 16 x 16, and 20 x 20). The case of the 
Hubbard model at U/t — 4 and L = 16 x 16 is also reported 
for comparison. Inset: Size scaling of the peak as a function 
of 1/L. 



long-range antiferromagnetic order disappears. In partic- 
ular, we find 5 C = 0.10 ± 0.01 and S c = 0.13 ± 0.02 for 
J/t = 0.2 and J/t — 0.4, respectively. 

At low doping, we have evidence that long-range order 
is always commensurate, with a (diverging) peak at X — 
(-7T, 7r) in the static spin structure factor, defined as 



) Q z Q z 



(13) 



This outcome is clear for all kinds of cluster considered, 
namely both for standard I x I and 45° tilted lattices. 
By contrast, close to the critical doping 5 C , we have the 
indication that some incommensurate peaks develop. Re- 
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markably, we do not find any strong doping dependence 
of the peak positions. We show the results of S(q) for 
the 16 x 16 cluster and J/t = 0.2 in Fig. [3l where the 
evolution of the peak as a function of the doping 5 is 
reported. By increasing the hole doping, the commen- 
surate peak at X reduces its intensity and eventually 
the maximum of S(q) shifts to a different k-point, i.e., 
(-7T, 7r — 2w/L). It should be stressed that this outcome is 
obtained only when the FN projection is applied to the 
lowest-energy ansatz contaning a sizable BCS parame- 
ter, and the FN calculation with a fully projected free- 
electron determinant cannot reproduce an incommensu- 
rate peak in S(q). Moreover, the variational wave func- 
tions always show commensurate correlations, see inset 
of Fig. [3l The strong dependence of this feature on the 
variational ansatz may also indicate that more accurate 
calculations are necessary to clarify this important as- 
pect of the phase diagram of the t—J model. In order to 
support the validity and the accuracy of our results, we 
have applied the same method to the Hubbard model at 
U/t = 4, where essentially exact calculations are avail- 
able for S l (q , )4i In this case, we have reproduced both the 
position and the intensity of the incommensurate peak on 
the 10 x 10 lattice. It is interesting to notice that, within 
the FN approximation, the intensity of the incommen- 
surate peak at U/t = 4 is much smaller than the corre- 
sponding one for the t—J model, see Fig. 2) This clearly 
indicates that in the t—J model the magnetic correlations 
are much more pronounced than the corresponding ones 
of the Hubbard model, possibly explaining the origin of 
the large extention of the antiferromagnetic region found 
in the t—J model. 

We now discuss whether these incommensurate spin 
correlations remain in the thermodynamic limit. For all 
the cluster sizes we considered, i.e., up to L — 20 x 20, 
the incommensurate peak in S(q) always appears at 
("7r, 7r — 2%/L), namely the closest k-point to X along the 
border of the Brillouin zone. This indicates that, in the 
thermodynamic limit, the peak should be located very 
close to X and it is not compatible with (n, 7r — 2ir5), 
found in cuprate materials^ Although size scaling ex- 
trapolations are not possible for a generic hole doping, we 
do not have evidence that the incommensurate peak di- 
verges, implying no incommensurate long-range order at 
finite doping concentrations. Nevertheless, once the com- 
mensurate magnetic order is melted, the ground state is 
characterized by short-range incommensurate spin corre- 
lations. In Fig. [U we show the results for J/t = 0.2 and 
S = 1/8, where different clusters with the same doping 
are available. Similar calculations with t' /t = —0.2 show 
the same qualitative behavior for the S(q). 

Coming back to the commensurate magnetic order 
close to half-filling, we stress that the pure t—J model 
shows robust antiferromagnetic correlations, with a crit- 
ical doping much larger than the one observed in the 
hole-doped cuprates materials, where the long-range or- 
der disappears at S c ~ 0.05.— This smaller value of the 
critical doping cannot be explained by reducing the an- 
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tiferromagnetic super-exchange J, given the fact that 
even for J/t = 0.2 we have that S c ~ 0.1. Besides dis- 
order effects that can be important in the undcrdoped 
regime^ one important ingredient to be considered in a 
microscopic model is the next-nearest-neighbor hopping, 
that was shown to have remarkable effects on both mag- 
netic and superconducting properties . 28 i 29 i 31 i 32 In partic- 
ular, in spite exact diagonalization calculations suggest a 
suppression of antiferromagnetic correlations for negative 
t'/i)22i more recent Monte Carlo simulations (also includ- 
ing a further third-neighbor hopping t") do not confirm 
these results, pointing instead toward a suppression of 
superconducting correlations 

In Fig. [5l we report the magnetization for J/t = 0.2 
and t'/t = -0.2. The first outcome is that the VMC 
results, though renormalized with respect to the case 
t' = 0, present a critical doping 6 C which is very simi- 
lar to the one found for the pure t—J model. By con- 
trast, the FN approach strongly suppresses the spin-spin 
correlations, even very close to half filling. In this case, 
the FN results have rather large size effects, that prevent 
us to extract a reliable estimate for the thermodynamic 
limit. However, it is clear that the antiferromagnetic re- 
gion is tiny and we can estimate that S c < 0.03. It should 
be emphasized that for t'/t = —0.2 the variational wave 
function is not as accurate as for the pure t—J model with 
t' = 0, but nevertheless the projection technique, even if 
approximate, is able to reduce the bias (e.g., the presence 
of a large magnetic order up to <5 ~ 0.1), showing the im- 
portance of alternative numerical methods to assess the 
actual accuracy of the simple variational approach. In- 
deed, we are confident that our FN results represent a 
good approximation of the true ground-state properties. 
On the contrary, the VMC calculations clearly show that 
the wave function @ overestimates the correct value of 
the magnetic moment. 
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C. Superconducting properties 

In the following, we want to address the problem of 
the superconducting properties of the Hamiltonian |TJ). 
In particular, we would like to obtain an accurate de- 
termination of the pair-pair correlations as a function of 
the hole doping and clarify the role of the next-nearest- 
neighbor hopping t' . The effect of such term has been re- 
cently considered by using different numerical techniques. 
DMRG calculations for rt-leg ladders (with n = 4 and 6) 
showed that the effect of a negative t' is to stabilize a 
metallic phase, without superconducting correlations. 29 
Moreover, improved variational Monte Carlo techniques 
suggested that t' could suppress pairing at low doping, 
whereas some increasing of superconducting correlations 
can be found in the optimal doping regime . 31 i 32 A further 
variational study^ suggested the possibility that a suffi- 
ciently large ratio t'/t can disfavor superconductivity and 
stabilize charge instabilities (stripes) near 1/8 doping. 

The pair-pair correlations are defined as 



0.025 



A">"(r) = S r ,A, 



(14) 



where SJ „ creates a singlet pair of electrons in the neigh- 
boring sites r and r + fi, namely 



-r,l S-+/i,r 



(15) 



For the first time, we implemented the forward-walking 
technique in order to compute true expectation values of 
the pairing correlations over the FN state: 



(A*"(r)> 



FN\ 



(16) 



Indeed, given the fact that A M,I/ (r) is a non-diagonal 
operator (in the basis of configurations defined before), 
within the FN approach the previous calculations^ were 
based upon the so-called mixed average, where, similarly 
to Eqs. (jS} and (j!2p. the state on the left is replaced by 
the variational one. Now, by using Eq. (fT5]) , it is possible 
to verify the fairness of the variational results against a 
much more accurate estimation of the exact correlation 
functions given by the FN approach. 

The superconducting off-diagonal long-range order im- 
plies a non-zero value of (A /i ' !y (r)) at large distance r. 
In the following, we consider the pair-pair correlation at 
the maximum distance and /i — v (parallel singlets) both 
for the VMC and the FN approximations and denote 
Pj = 41im r ^ 00 (A y ' y (r)). It is worth noting that, as far 
as the superconducting correlations are concerned, there 
is no appreciable difference between the results obtained 
with and without the antiferromagnetic order parame- 
ter and the long-range spin Jastrow factor. The results 
for the pure t—J model are reported in Fig [6) where we 
report two different values of the antiferromagnetic cou- 
pling, i.e., J/t = 0.2 and J/t = 0.4. In this case, VMC 
and FN calculations arc in fairly good agreement, giv- 
ing a similar superconducting phase diagram. Interest- 
ingly, the optimal doping, i.e., the doping at which the 
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maximum in the pair-pair correlations takes place, oc- 
curs in both cases at S ~ 0.2, whereas the actual value of 
the correlations is proportional to J/t. At high doping, 
where antiferromagnetic fluctuations play a minor role, 
the behavior of the pairing is unchanged when J is var- 
ied. Although in this region there are some size effects, 
we can safely estimate that superconductivity disappears 
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around 6 ~ 0.35 and S ~ 0.4 for J/t = 0.2 and J/t = 0.4, 
respectively. 

It is worth noting that the density Jastrow term (j4|) 
is very important to obtain an accurate estimation of 
the pairing correlations. Indeed, the variational results, 
based on the simple wave function with BCS pairing and 
the on-site Gutzwiller projector (but without the long- 
range Jastrow factors) overestimate the pairing correla- 
tions at optimal doping by at least a factor two. Remark- 
ably, the FN approach is able to correct this bias, provid- 
ing approximately the same results as the one obtained 
starting from the wave function with the long-range Jas- 
trow factor, see Fig. [6] This fact shows that the FN 
method is particularly reliable for estimating the pairing 
correlations. 

The inclusion of the next-nearest-neighbor hopping in- 
duces sizable modifications in the pairing correlations, 
though the qualitative dome-like behavior remains un- 
changed, see Fig. [JJ At low doping there is a sizable 
suppression of the superconducting pairing, particularly 
evident after the FN projection, see Fig. [8j Indeed, while 
for the pure t—J model we clearly obtain a linear be- 
havior of the pair-pair correlations with indicating 
a superconducting phase as soon as the Mott insulator 
is doped, in the case of a finite t' , the FN results could 
be compatible with a finite critical doping, below which 
the system is not superconducting. This outcome is in 
agreement with earlier Monte Carlo calculations done by 
one of us^i where it was suggested that the extended 
t—J model with hoppings and super-exchange interac- 
tions derived from structural data of the La2Cu04 com- 
pound could explain the main experimental features of 
high-temperature superconducting materials, with a fi- 
nite critical doping for the onset of electron pairing. 

Remarkably, from S ~ 0.1 to S ~ 0.4 there are huge 
size effects. Though, for S ~ 0.3, small clusters, e.g., 
L = 98, indicate stronger pairing correlations than the 
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pure t—J model without t', larger clusters point out a 
large reduction of Pj. Nonetheless, we have a rather 
clear evidence that for S ~ 0.3 there is a finite super- 
conducting order parameter in the thermodynamic limit, 
see Fig. |H1 This strong reduction of the superconducting 
correlations is a very interesting effect, demonstrating 
that the superconducting wave function (even if supple- 
mented by magnetic order) deteriorates its accuracy by 
increasing the value of t', that could eventually stabilize 
competing phases with modulation in the charge distri- 
bution and/or a magnetic flux through the plaquettes4^ 
However, for t'/t = —0.2, the homogenous variational 
ansatz ^ provides a lower energy when compared to 
the one used in Ref. I45l4£ At present, the most accurate 
FN calculations based on the lowest-energy variational 
ansatz, do not show any tendency towards charge in- 
homogeneities for S < 0.4. This outcome is important 
because in principle the FN approach can spontaneously 
induce charge-density wave modulations in the ground 
state, even when the variational wave function before the 
FN projection is translationally invariant. 12 



IV. CONCLUSION 

In this paper, we considered the magnetic and su- 
perconducting properties of the t—t' — J model within 
VMC and FN approaches. We showed that for t' = 
the ground-state properties can be accurately reproduced 
by a state containing both electron pairing and suit- 
able magnetic correlations, namely a magnetic order pa- 
rameter in the mean-field Hamiltonian that defines the 
fermionic determinant and a spin Jastrow factor for de- 
scribing the spin fluctuations. In this case, we obtain a 
rather large magnetic phase, with a critical doping that 
slightly depends upon the super-exchange coupling J, 
i.e., S c = 0.10 ± 0.01 and S c = 0.13 ± 0.02 for J/t = 0.2 
and J/t = 0.4, respectively. The superconducting corre- 
lations show a dome-like behavior and vanish when the 
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Mott insulator at half filling is approached. Interestingly, 
compared to the RMFT that predicts a quadratic behav- 
ior of the pair-pair correlations as a function of the doping 
5, here we found that a linear behavior is more plausible. 

Then, we also reported important modifications due to 
the presence of a finite ratio t'/t. The first effect of this 
further hopping term is to strongly suppress antiferro- 
magnetic correlations at low doping, shifting the critical 
doping to 0.03 for t'/t — —0.2. This is a genuine effect 
of the FN method, since, within the pure variational ap- 
proach, though the spin-spin correlations are suppressed 
with respect to the case of t' = 0, the values of the crit- 
ical doping for these two cases are very similar. Most 
importantly, the presence of a finite value of the next- 
nearest-neighbor hopping has dramatic effects on the su- 
perconducting properties. At small doping, i.e., 5 < 0.1 
there is a sizable suppression of the electronic pairing, 
possibly pointing toward a metallic phase in the slightly 
doped regime, as previously suggested by using improved 
Monte Carlo techniques. 44 Moreover, for 0.1 < S < 0.4, 
though small lattices seem to indicate an increasing of su- 
perconductivity compared to the pure t—J model, larger 
clusters show huge size effects that strongly renormalize 
the pairing correlations at large distance. However, for 
the value of t' considered in this work, we are rather con- 



fident that superconducting off-diagonal long-range order 
takes place in a considerable hole region. In any case, the 
huge renormalization of the electronic pairing for 6 ~ 0.3, 
together with the fact that the FN results are very differ- 
ent from the VMC ones based on a wave function contain- 
ing pairing (and magnetic order at low doping), suggests 
the possibility of the existence of a non-superconducting 
phase (with magnetic fluxes and/or charge order) that 
could be eventually stabilized by further increasing the 
ratio t'/t. 

In conclusion, the main qualitative features of the 
cuprate phase diagram appear rather well reproduced by 
the t—t'—J model with a sizable next-nearest- neighbor 
hopping. However, we do not find a sizable enhancement 
of the pairing correlations by increasing the ratio t'/t, 
that looks in contradiction with the empirical relation 
between t'/t and the value of T c .— However, we have to 
remark that we only considered ground-state properties 
and we cannot evaluate T c , whose relation with the pair- 
ing correlations may be highly non trivial, especially in 
a strongly-correlated system, like the t—J model. This 
issue certainly deserves further work. 

We thank useful discussions with S. Bieri, M. Capello, 
and D. Poilblanc. This research has been supported by 
PRIN-COFIN 2005 and CNR-INFM. 
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